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SUMMARY 



This report presents a survey of previous theoretical 
and experimental investigations on wing flutter covering 
thirteen cases of flutter observed on airplanes. The di- 
rect cause of flutter is, in the majority of cases, at- 
tributable to (mass-) unbalanced ailerons. 

Under the conservative assumption that the flutter 
with the phase angle most favorable for excitation occurs 
only in two degrees of freedom, the lowest critical speed 
can be estimated from the data obtained on the oscillation 
bench. Corrective measures for increasing the critical 
speed and for definite avoidance of wing flutter, are dis- 
cussed. 

I. INTRODUCTION 



The forced oscillations on airplane wings are oscil- 
lations created solely by the air stream and have as a 
rule nothing to do with the vibrations set up by the in- 
ertia forces of the engine. They are therefore best des- 
ignated by the term 11 flutter" since they revert to the 
same underlying causes as the fluttering of a flag. 

Flutter starts at the so-called "critical speed, 11 
which depends chiefly on the oscillation frequency and on 
the wing chord. The lower the frequency and the smaller 
the chord, the lower the critical speed will be. The os- 
cillation frequency of a wing, in turn, depends on the 
stiffness and on the mass of the wing. 

Flutter in an air stream is possible only when a 
plate - in whole or in part - is free to rotate about at 



" Augenblicklicher Entwiclrlungsstand der Frage des Flugel- 
flatterns." Luf t f ahrtf or schung , October 3, 1935, pp. 
193-209. 
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least two axes or, which is the same, has at least two de- 
grees of freedom of oscillation. 

A wind vane of sheet metal made to rotate about one 
axis only does not flutter. However, when the flagpole is 
not rigid hut free to swing laterally, thereby pivoting the 
vane about an axis "below it and parallel to the wind, flut- 
ter is possible. If the vane is of cloth rather than metal, 
it can turn about infinitely many axes and is therefore 
particularly susceptible to flutter. So also is a wind 
vane made of two pieces of sheet metal hinged together, be- 
cause then the flagpole and the hinge line between the two 
pieces form the two axes of rotation,. 

A similar condition exists when mounting a rudder E 
with tab H to a practically rigid fin J (fig- l) . The 
two axes of rotation are A x and A s . With such an ar- 
rangement flutter has actually been observed (reference 1). 

Far more importance, from the practical point of view, 
attaches to the case of an airplane wing fitted with an 
aileron. When oscillating, the wing, turns about some nodal 
axis which may, for instance, coincide with the wing center 
line or the axis of the strut connections* Besides, the 
aileron itself can turn about its hinge. 

The first records of wing flutter go back to the ear- 
ly days of flying, when the lateral control obtained by 
twisting the wing tips, was abandoned in favor of the aile- 
ron-control method 

During the World War several cases Occurred where 
flutter caused the ailerons to break and tear off. Like- 
wise, almost all cases observed later on disclosed upon 
investigation, that the ailerons were the cause of the ac- 
cident. Even a rigid plate can flutter, as stated above, 
when free to rotate about two axes, If the wing tip bends 
and twists simultaneously, it can flutter even without ai- 
lerons, although this case is much less frequent than the 
one described first. 

In the following, the results of past investigations 
on wing flutter are given without resorting to mathematical 
deductions, while one section contains a discussion of the 
theoret ical relat ions. 
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II. DEVELOPMENT OF METHODS OP ANALYSIS 



The exploration of the causes of wing flutter is 
marked "by the diversity of methods employed with a view to 
obtaining technically useful solutions of this extremely 
complicated problem. 

1. Theorem of Linear Differential Equations 

for Steady Aerodynamic Forces 

In the first flutter investigations, the air loads on 
the oscillating wing were assumed to "be steady and depend- 
ent on the dynamic angle of attack, the dynamic angle of 
attack "being defined as the angle "between the wing chord 
and the momentary direction of motion of the oscillating 
wing. Some authors also took into account the lift due to 
dynamic profile camber. ' A wing oscillating about some ax- 
is, while its wing chord describes a curved surface line 
in flight, is identical with a wing in steady flight whose 
profile curvature changes at measured intervals. 

This substitution is, in fact, strictly correct. 
Even these elementary assumptions afford a physical expla*- 
nation of the phenomenon of flutter through a system of * ' 
linear differential equations, the number of which depends 
on the number of degrees of freedom. Flutter is possible 
whenever undamped oscillations of constant amplitude, i.e., 
harmonic oscillations, are possible. Routh' s discriminant 
thereby served as a criterion from which the critical speed 
may be computed. . 

The first calculations of this kind were made by 
Blasius in June 1918, at the request of the Inspection 
Section of the German Air Corps (reference 2), incident to 
the investigation of the flutter on the lower wing of the 
Albatros D3 biplane which, having only one spar, was of 
low torsional stiffness. There were no ailerons on the 
lower wing. The accidental circumstance which prompted 
the investigation of that particular case at all, was due 
to the fact that at that time the significance of the ai- 
leron as promoter of flutter, was not sufficiently appre- 
ciated. Similar investigations were subsequently made by 
v, Baumhauer and Zoning, Bairstow, Frazer and Duncan, 
Blenk and Liebers, Hesselbach, and were extended to include 
oscillating ailerons (references 3 to 14), 
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Indeed, the calculation of the simple" elastic-mass os- 
cillations of an airplane wing on the oscillating "bench, 
stipulated a number of simplifying assumptions. Other sim- 
plifying assumptions consisted in disregarding the energy- 
consuming, unsteady system of vortices and the premises of 
material damping proportional to the rate of deformation. 
The accuracy of such calculations therefore is, as a rule, 
quite small. numerical agreement between calculation and 
experiment has "been obtained only in cases where the as- 
sumptions could be made to fit the particular case. 

Agreement was more readily obtainable on cantilever 
than on braced wings. At first it was believed that can- 
tilever monoplanes were particularly susceptible to flutter, 
but subsequent experience proved otherwise. 

One important result was the following rule: The mass 
axis of the wing shall lie ahead of the classic axis if 
feasible; the aileron e.g. shall lie in its hinge axis in 
order to avoid flutter. 

2. Calculation of Vortex Separation 

Whereas in the early stages of development, wing flut- 
ter was treated as a mechanical problem, the aerodynamical 
side now received more attention and it was attempted to 
trace the source of the unsteady lift of the oscillating 
wing and the correlated separation of vortices, at least 
for the case of two-dimensional flow. 

The problem of the oscillating wing was first attacked 
by Birnbaum (references 15 and 16). He introduced the im- 
portant concept of the reduced frequency CO, which is tt 
times the ratio of wing chord to wave lengths. If n is 
the oscillation frequency (in minutes), t, the wing chord 
(in meters), and v, the flying speed (in kilometers per 
hour), the reduced frequency is: 

0) = 0.06 tt S£ (1) 

The air loads on the oscillating wing are functions of this 
nondimensional parameter. Following the example of Prandtl, 
Birnbaum replaced the wing by a system of bound vortices 
and postulated that the sum of bound and free vortices must 
remain constant with time; he obtained an equation which he 
could solve for small values of the reduced frequency 
0) < 0.12. Beyond this point his development was not con- 
vergent . 
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This solution, however, was -unsuitable for the eluci- 
dation of the problem under consideration, "because wing 
flutter always occurs at materially higher values of oo . 
Kussner found the general solution of Birnbaum 1 s equation 
and extended it to include the case of the co-oscillating 
aileron (reference 17). Since the creation of harmonic 
oscillations is always considered as the oscillation cri- 
terion, it seemed natural to write the equations from the 
very first for harmonic o sc illat ions which offers the 
added advantage of utilizing the labor-saving method of 
complex presentation. The oscillation criterion then is 
the disappearance of the complex denominator determinant, 
which yields two equations for calculating the oscillation 
frequency and the critical speed. This obviates the use 
of the linear differential equation and Routh f 3 discrimi- 
nant* One particular advantage accruing from the use of 
the harmonic oscillation is that the material damping can 
be introduced in a simple and physically correct manner as 
phase difference of the elastic force. This possibility 
does not exist with the linear differential equation, where 
it is even necessary to make a physically incorrect assump- 
tion of the damping in order to obtain a linear equation. 



On this basis it was then possible to calculate sever- 
al examples of an oscillating flat plate in order to elu- 
cidate systematically the influence of mass distribution, 
elastic forces, and material damping. It was found that 
with two degrees of freedom -r bending and torsion - the 
critical speed depends chiefly on: 

1. The torsional oscillation frequency of the wing. 

2. The backward position of the e.g. of the wing. 

3. The material damping. 

The result of material damping is that flutter is 
possible only up to certain maximum In oscillations 

at higher w, the energy obtainable from the air stream 
would become inadequate for compensating the damping losses. 
This rule holds not only for the two degrees of freedom un- 
der discussion.- bending and torsion - but is of general 
validity, as will be shown later. 

Theoretically the effect of material damping is so 
much greater, as the ratio of bending stiffness to tor- 
sional stiffness is higher, which is approximately equiv- 
alent to the ratio of wing chord to length of overhang. 
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As the total damping is not accurately determinable, no 
fairly close agreement can "be expected between calcula- 
tion and observation except for cantilever wings of high 
aspect rat io . 

Similar investigations have been made in England, 
Glauert calculated the unsteady air loads on an oscillat- 
ing wing for the two degrees of freedom - bending and tor- 
sion (reference 18). He proceeded from H. Wagner's con- 
cept of the area of discontinuity; but as his numerical 
calculations e;ctend only to (0 = 0.5, they are insuffi- 
cient for the mathematical treatment of flutter, with its 
much higher to values. 

Duncan and Collar extended the calculation to a wing 
oscillating with increasing amplitude (reference 19). 
Lately, Theodorsen has calculated the air forces on an os- 
cillating airfoil (reference 19a). 

3. Model Experiments 

Model experiments are another means of investigating 
wing flutter, but if such model tests are to afford prac- 
tical conclusions the models must be constructed dynamic- 
ally similar. Dynamic similarity is the more difficult to 
attain as the model scale, i.e., the model, is smaller. 
Since the model scale depends moreover on the jet diameter 
of the available wind tunnel, the dynamic similarity was 
disregarded at first and the simply constructed model 
wings were mounted in the air stream to a wall represent- 
ing the plane of symmetry of the wing (references 9, 22, 
and 23). Such models were sufficient for exploring the 
effect of e.g. position, damping, and mass unbalance of 
the aileron. But the values of the reduced frequency ob- 
tained in these tests are considerably less than the ex- 
perimental values cited below. 

The British have investigated a great number of ac- 
tual cases of flutter besides model testing since 1925, 
and have shown great skill in their choice of assumptions 
which afforded agreement between calculation and observa- 
tion (references 7 to 13). Model experiments were fre- 
quently used as basis for computing the still unknown 
damping forces, the linear differential equations forming 
the starting point, while South 1 s discriminant was ex- 
pressed as determinant, whereby some fields of the deter- 
minant remained empty. 
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The method of calculation given in reference 17 was 
checked at the D.V.L. by wind-tunnel experiments on model 
cantilever wings, which were, of course,, fairly heavy and 
as a result, oscillated at a lower reduced frequency 

0,3 (reference 22). The observed critical speed on 
five model wings was from 14 to 24 percent higher than the 
theoretical, which may be attributed to the flow being 
other than two-dimensional and to the energy absorption of 
the disregarded trailing vortices. 

Subsequently two dynamically similar models of the 
He 60 type were constructed ,at 1:5.6 scale with a span of 
2,4 meters. These model tests were intended to trace the 
cause of the accident described elsewhere and to test the 
efficacy of certain structural changes with a view to pre- 
venting flutter. The problem was solved, although a num- 
ber of unexpected difficulties were encountered in this 
first attempt at constructing dynamically similar models. 
The highest reduced frequency obtained in the tests was 
0.76, a figure which is fairly close to the probable co = 
0.93 at the time of .the accident* Since complete dynamic 
similarity is not attainable and the model usually has more 
damping than the full-scale wing, the expected co value 
for the model will in any case be less than for the full- 
size wing. 

4. Statistical Investigation 

Admittedly, the methods of investigation described so 
far suffice to explain observed cases of flutter and to 
prove the underlying causes of such flutter, wherein the 
actually observed critical speed always constituted a 
check on the correctness of the assumptions. But these 
methods did not lend themselves to computing the critical 
speed on a new type of airplane within a fair degree of 
accuracy, particularly when applied to braced wings. Any 
further analytical treatment of flutter was precluded, 
since it was impossible to compute the purely elastiq os- 
cillations of a wing on the stand with a reasonable amomt 
of paper work, unless the construction was fairly simple, 
such as monospar, cantilever wings. As a result it was at- 
tempted to establish a simple dimension rule, suitable for 
practical use by the designer, to prevent wing flutter due 
to torsion within the normal speed range. 

Since the wing mass and its backward e.g. position 
are little amenable to influence, aside from the fact that 
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the naterial damping should also be considered as predeter- 
mined, the only valid means for raising the critical speed 
is the torsional stiffness. The German design specif ica- . 
tions carried a provision for torsional stiffness as far 
back as 1918 for army airplanes. The angle of twist in 
the terminal dive was not to exceed 5°; this was reduced 
to 3.5° in the 1926 design specifications. However, this 
was primarily with a view to static torsional stability of 
the wing rather than to wing flutter. With the increasing 
use of airfoils with fixed c.p. , this requirement became 
useless. 

The 1930 specifications contained a rule of thumb for 
torsional stiffness, based on a few theoretical examples 
and similarity considerations (reference 17).* 

r (y) = % ft ^ * Po v k s F (yf ( k ^ S > < 2 > 

In this formula, k was at first put at k = 0.12 to 
0.24, but subsequent calculations brought about a change 
to k = 0.5 (1934 design specifications). With this as- 
sumption it is already very probable that the true critical 
speed lies above that given in formula (2). It was there- 
fore permissible to introduce the terminal diving speed 

v in formula (2) . 
c 

It is worthy of note that this formula, originally 
merely intended for the degrees of freedom - wing torsion 
and bending - proved practical also for a number of air- 
foil-aileron combinations, because the observed maximum 
values of the reduced frequency for this type of oscilla- 
tion are of the same order of magnitude as the frequencies 
stipulated for wing bending and torsion. 

Hoxbee Cox checked formula (2) against ten actual 
flutter cases (references 20 and 2l). He applied torque 
U& at the wing tip A- A , measured the angle of twist cp , 
and computed therefrom the constant 

k , = JJd, lSL (3) 

p o k (y) 

Its numerical values are given in tatle I. 

*Pt>r symbols, see section IV, 1. 
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TABLE I. Torsional Stiffness and Flutter 



Type 


m 3 


.*« 
in 


V k 
kn/h 


k' 


k 


Gloster " Gamec o ck" 


6. 68 


1 . 60 


258 


0 .135 


0 . 40 5 


G In q t " Gn mer o fk" in 
pull-out 


S. 68 


1. 60 


403 


.074 


.224 


Gloster " Sor co ck" , wood 


6.58 


1.60 


290 


.10 6 


.318 


Gloster "Gorcock" , metal 


6.58 


1.60 


217 


.257 


.771 


Short "Satellite" 


7 .87 


1. 68 


145 


.033 


.066 


Gloster "Grebe" 


6.54 


1, 60 


2 58 


.116 


.348 


Desoutter Mark II 


6.97 


1.55 


22 5 


.295 


. 590 


Mar t insyde F* 4 


6.08 


1.68 


323 


.098 


.294 


DeHavilland "Puss Moth" 


8.80 


. 1.83 


314 


.129 


.258 


Simmond "Spartan" 


5.02 


1.37 


274 


.072 


.216 



m 2 X 10.7639 = sq.ft. kn/h X .62137 = ni./hr. 



The characterization of the torsional stiffness solely 
through angle of twist at the tips is a rather summary pro- 
cedure. Consequently, the k 1 values scatter considerably. 
If the increase in angle of twist at the tip of a monoplane 
wing is twice as great n,s the mean value over the whole 
wing and three times as great for a biplane wing, then the 
k values given in the last column of table I are compara- 
ble to the mean value k = 0.35 (formula (3)). Only two 
values lie above the maximum value of 0#5 stipulated in 
the 1934 design specifications. It seems reasonable to 
assume that these two cases at least involve flutter with 
wing flexure and aileron motion. Unfortunately the Brit** 
ish report fails to give the modes of oscillation and the 
flutter frequencies. Index values for the bending stiff- 
nesses were established in a similar manner. However, it 
serves no useful purpose to analyze these makeshift dimen- 
sion rules, because section IV contains a method which af- 
fords a better estimate of the critical speed. 



10 



U.A.C.A, Technical Memorandum Ho* 782 



The stiffness formulas are makeshift substitutes for 
the calculation of the purely elastic oscillation frequen- 
cies of a cellule - a calculation which is often quite 
difficult. This difficulty may he overcome by subjecting 
the finished airplane to a static oscillation test. To 
this end the, airplane is elast ically mounted, an unbalanced 
rotor is attached "below the fuselage and driven at varying 
speeds by an electric motor, the mode of oscillation and 
the frequency "being recorded iji resonance conditions. 

Even when the data of such oscillation tests are 
available, it is still extremely difficult and tedious to 
analyze the critical flutter speed for the three degrees 
of freedom - bending, torsion, and aileron motion - because 
the calculation still contains important simplifying as- 
sumptions, especially that of two-dimensional flow, as a 
result of which the possibility of error should not be un- 
derestimated, It is true, however, that this error is 
usually on the safe side, as shown by the comparison be- 
tween calculations and model tests mentioned above, be- 
cause any damping, neglected in the calculation, will 
raise the critical speed. In such a calculation, made 
with the utmost care, for the braced He 9a monoplane a re- 
duced frequency of 0) = 1,13 was established, the possible 
error being estimated at -20 percent. The chief drawback 
of the operation lies in the physically correct terras for 
the complex determinant rather than in the evaluation of 
the determinant ■ 

Presumably no substantially greater mathematical accu- 
racy can be obtained even after the calculation has been 
improved and refined, because flutter does not always start 
at the sane speed even in the wind-tunnel test. The turbu- 
lence of the air stream, the angle of attack of the wing, 
and accidental snail differences in the hinge friction - 
all have some influence. Past experience has been that 
flutter often starts in gusty weather, from which it may 
be concluded that gust shocks have overcome the initially 
excessive friction forces. 

Once flutter has started - in this or some other man- 
ner - it frequently continues until the pilot has reduced 
the speed to two thirds or less of its original value. 
Possible causes, for this are: rupture of the aileron con- 
trol cables, the consistently smaller proportion of hinge 
friction to the total damping as the amplitude increases, 
and lastly, the effect of change in angle of attack. 
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- Another fact should mentioned in this connection. 
At very low ampl itudes the laws of potential flow do not 
hold "because then the viscosity of the air is no longer 
negligible. Consequently, the air forces are smaller for 
very small amplitudes than they should "be according to 
the potential theory,' and therefore do not induce flutter. 
This effect was observed by Birnbaum (reference 15, p. 
292), It apparently is a boundary-layer effect. The wing 
flops around, so to speak, in its boundary layer, without 
encountering any resistance. 

In flight free from oscillations and at uniform speed 
through still air, a wing could exceed its critical speed 

by any. amount without starting to flutter. It would take 
a shock of a certain minimum size, e.g., a gust shock, to 
start flutter which, on the other hand would, of course, then 
be extremely violent. ilotable in this connection is the 
fact that, flutter has often been observed during or direct- 
ly following a pull-out from a steep dive, particularly 
in vicious case's.' Iri 'a. normal , mild p'uir-oAit from high 
speeds, only small changes of angle of attack are po-ssible. 
It is improbable that the' quotient d c a /da, on which the 

air forces depend, changes very materially within such a 
small range of angle of attack. One may suspect, there- 
fore, that the disturbance of the. boundary layer during 
transition from gliding to pull-out or pull-out to level- 
off was the trigger' effect in these, cases. 

Summing up these facts deduced from experience and 
considering in particular the great amount of time' re- 
quired for the calculation, which is not justified by the 
small degree of accuracy, one comes to the conclusion that 
the analytical method, while adequate for explaining the 
fundamental relations, is scarcely suitable for the predic- 
tion of the critical speed of a new type of airplane. 

Once a physical process i s no longer amenable to ana- 
lytical treatment because it contains variables which can- 
not be observed and numerically defined, then it must be 
explained statist icalljr , ' based on. a large number of obser*- 
vations. This statistical method, indicated during the 
formulation of the stiffness formula (2) , can now be ap r 
plied in a more comprehensive manner to the problem of 
wing flutter, because within the last few years a number 
Jf. cases' of flutter have been ihvo st igated - in detail , even 
though, this number is as yet not very large from the point 
of view of statistical research 
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The most important parameter introduced in the' analyt- 
ical treatment is the reduced frequency. An attempt was 
therefore made to determine the reduced frequency in the 
observed cases of flutter. 

One may differentiate "between "mild" and "vicious" 
cases. In mild cases flutter occurs with small amplitudes 
which are well "below the ultimate strength of the wing. 
The flutter usually stops at a speed slightly "below that 
at which it started, so that the flutter may "be stopped 
very quickly "by pulling the stick "back. These mild cases, 
while few in number, can "be demonstrated with comparative- 
ly little danger and are therefore suitable for flutter 
investigations in free flight. A test of this kind made 
on the He 45c , is described elsewhere in the report. The 
recorded air speed, frequency, and mode of oscillation in 
flight affords the true value of the reduced frequency and 
the ratio of the amplitudes for each degree of freedom*' 

This determination is more difficult in the vicious 
cases. In these cases flutter is, in a way, actually de- 
layed by the very causes cited above and does not start 
until the theoretical critical speed has been exceeded; 
then, however, it begins with such violence as to cause 
failure of the wings or ailerons. If the airplane is still 
able to land, it is repaired after the flight and subject- 
ed to an oscillation test. The dangerous mode of oscilla- 
tion is that at which the lowest frequency is accompanied 
by torsional oscillations of the wing or aileron for the 
reason that, aside from wing flexure, it requires one of 
these two degrees of freedom to give increasing amplitudes. 
However, this does not imply that flutter must occur at 
the frequency observed in the oscillation test, because the 
air forces existing during flutter may modify the mode of 
oscillation and the frequency. In particular, a differ- 
ence in phase angle is always to be expected between bend- 
ing and torsion, because it is only under these conditions 
that the energy for increasing the amplitudes can be taken 
out of the air stream, Even so, the oscillation test af- 
fords a certain basis, which is the more reliable as the 
resonance condition appearing in the oscillation test is 
more definitely expressed; i.e.* as the damping is small- 
er. (See section IV, 4.) 

If the airplane is destroyed by. th.e accident, another 
airplane of the same type will be subjected to the oscilla- 
tion test. The flight speed at the time of the accident 
can rarely be given very accurately for obvious reasons. 
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In serious accidents one nay have to rely on statements of 
ey ewit nesses on the ground. 

The possibility of errors introduced when determining 
the reduced frequency is therefore great for the vicious 
cases. It is also necessary to decide whether the calcula- 
tion of the reduced frequency is to he effected at the 
spe ed at which f lut ter started or at whi ch it stopped . 
But since the start of flutter is decisive for flight op- 
eration, the speed at incipient flutter is customarily pre- 
ferred* In this manner the observed values discussed in 
the next section have been obtained. Disregarding the 
possible errors, they range between oo = 0.58 and od = 1*14, 
from which it appears that the reduced frequency in new 
types of airplanes will not exceed 0) h = 1.14. 

Testing an airplane on the oscillation bench and ob- 
serving the dangerous mode of oscillation with the fre- 
quency n, the. lowest possible value of the critical speed 
can be roughly estimated on the basis of the assumed maxi* 
mum value of the reduced frequency. If t m is the 

mean chord of the outer part of the oscillating wing, the 
lowest possible value of the critical speed is 



Obviously such a statistical .appraisal is wo rthless ' unl ess 
the particular type of airplane is not substantially dif- 
ferent from all the airplanes which showed flutter in the 
indicated range of reduced frequencies by having incorpo- 
rated special features which minimize flutter hazard. 
When these investigations on flutter were started, the 
probability of finding such a type of airplane was very 
small, but in t ime there -wi 11 he an ever- increasing number 
of types on which such preventative measures may be effect- 
ed with at least the partial success of lower reduced fre- 
quency. This being so, the rough statistical estimate may 
be replaced by an. improved method (section IV) which per- 
mits the inclusion of proved preventative measures. 

For mass-balanced ailerons or wings without ailerons, 
the lowest possible critical speed is higher, and the re- 
duced frequency consequently lower, than the maximum value 
given above. Practical data are very scarce on this sub- 
ject, because inall cases of flutter described hereinaf- 
ter, the ailerons contributed to the growth of oscillations; 




m 



km/h 



(4) 
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at least, it was impossible to state whether in one case 
or the other, flutter would also have occurred if the aile- 
rons had "been rigidly connected to the wings. The one case 
of flutter without aileron which had "been definitely estab- 
lished, prompted the first investigation of flutter without 
aileron (section II, l) . 

llaturally, as the speed of airplanes increase, the re- 
gion of flutter without aileron will also he reached more 
frequently, and any appraisal of the critical speed based 
on the static oscillation must allow for this possible mode 
of oscillation also. 

III. RESULTS 
1. Analysis of Observed Cases 

a ) Braced D P 9 (references 4 and 54). - The s t rut i s 

short, so that a long overhang exists. This model devel- 
oped two cases of , vicious flutter in the spring and autumn 
of 1925, starting during pull-up from a steep glide at 
about 180 km/h# In one case it led to complete fracture 
of the wing; in the other, to fracture of the ribs in the 
overhang and of the aileron control cables. In gusty „ 
weather it started a slight flutter at 135 km/h. 

After the wing was mounted on a rigid test frame, it 
showed a flexural oscillation frequency of 548/min. , and a 
torsional oscillation frequency of 494/min. , The frequency 
of the free oscillation may be rated at 520/min. The wing 
chord was 1.5 m; the aileron chord 0.32 m; and the aileron 
eg. was 126 nm behins the hinge line. The reduced fre- 
quency is 

co = 0.1885 =0.82 
180 

b) Braced He 8a no nop 1 an e . - This airplane crashed in 

the fall of 1928, due to fracture of the wings during an 
exhibition flight. From the reports of eyewitnesses, it 
seems quite safe to conclude that flutter was the cause. 
The flight speed is estimated at 350 km/h. An airplane of 
the same type was tested on the oscillation stand. The 
dangerous mode of oscillation lies probably at 540/min., 
and has a nodal line running from the rear strut fitting 
toward the point where the curved tip joins the straight 
leading edge (fig. 3). 

m X 39.37 = in. on X .03937 = in. kn/h X .62137 = mi./hr. 
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The wing chord was 3.0 m, the aileron chord 0,3 in, 
the aileron eg; was 50 mm "behind the hinge line. The re- 
duced frequency* is 

u> = 0.1885 = 0.87 

3 50 

c) Bra ce d _L _ 7_8 Jb ipl_ an c . - Th i s mo d e 1 , of which quite a 

number had "been built, had often been dived at 350 km/h,' 
when in I/lay 1930, it developed a case of nild flutter while 
flying at a speed of about 210 km/h. It started with an 
oscillation of the strut between the lower and upper aile- 
rons at great amplitude, then the wings fluttered so se- 
verely that the pilot was 4 unable to hold the stick. As 
soon as the pilot cut his 4 speed, the oscillations died 
out. The dangerous mode lies at 860/min, The lower wing 
oscillates in bending, the nodal line being near the strut 
fittings. 

The aileron connecting strut shows severe lateral de- 
flections which cause the upper-wing ailerons to oscillctte 
in torsion ( fig. 4) . 

The mean chord of the overhang of the lower wing is 
1.36 m; the aileron chord from hinge line to trailing edge 
is 115 mm; the aileron e.g. is 23 mm behind the hinge line. 
The reduced frequency is 

co = 0.1885 -~2fo~ ~ = 1.05 

d) Unbraced. ,_He 60 biplane.- This model is a rather 

less conventional design. The lower wing is braced against 
the floats while the two struts on each side reach only to 
the front spar. Ho wire bracing is used between the wings. 
While in other versions of this type the spars had been 
made of wood, this particular type (He 60) utilized steel, 
providing the same strength for the same spar height. The 
ratio of Young* s modulus to ultimate strength for steel 
being substantially lower than for wood, it assured low 
natural frequencies of the wing, In addition, the aileron 
system had an unbalance of 75 cm kg. Apprehensions were 
therefore voiced from the very beginning that flutter might 
occur at speeds lower than the prescribed diving speed of 

3 65 km/h. 

In the attempt to reach the prescribed diving speed, 
the airplane crashed in December 1931, as a result of a 
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torn upper wing. According to the testimony, of eyewit- 
nesses, it was a case of dangerous flutter. The speed was 
estimated at 350 kru/h. 

The dangerous node lies at 780/nIn. The nodal line 
of the upper wing was in the overhang close to and almost 
parallel with the leading edge in the inner "bay; at approx- 
imate wing center it runs parallel to the wing axis (fig* 
5). The wing chord is 2*2 m; the aileron chord 0.4 m. 
The reduced frequency is 

0) = 0.1885 = 0.93 

o) Braced He 46c biplane,- The He 4Sc is a braced "bi- 
plane with a snail lower wing developed from a high-wing 
monoplane. After a long period of service, it finally re- 
vealed a mild case of flutter at 260 kn/li which, however, 
disappeared immediately as soon as the speed was reduced. 
It was therefore decided to obtain some oscillation photo- 
graphs in flight with this airplane, taking, of course, 
proper precautions* The records showed the flexural oscil- 
lation of the lower wing, coupled with turning of the un- 
balanced aileron system as the cause of the increase; the 
aileron c#g* was 52 mm behind the hinge line. The aileron 
chord of the lower wing is 335 mm, and that of the upper 
wing, 500 nn. ; the wing chord is 1.4 m on the lower, and 
2.0 m on the upper wing. 

The oscillation test disclosed between 520 and 755/ 
nin. > a series of antisyrnmetrical oscillation modes of the 
whole cellule about the longitudinal, verticals-transverse, 
axes, accompanied in part by severe aileron motions (aile- 
ron control by means of torque tubes) • The remarkable fea- 
ture is that these modes do not induce flutter. This may • be 
attributable to a slight mass- coup ling with the aileron os- 
cillation as a result of the small amplitude of the aile- 
ron hinge line and the shifting of the location of the 
nodal line closer to the trailing edge. Possibly the damp- 
ing of the antisyrnmetrical oscillations of the whole cel- 
lule is greater. The first symmetrical natural bending 
frequencies of the wings lie between 815 and 895/min. The 
nodal line of the lower wing lies at 50 percent or more of 
wing chord forward of the leading edge (fig. 6). The reso- 
nance conditions are not pronounced. 

But the flutter frequencies recorded with the opti- 
graph, lie in this range. Flutter started with a frequency 
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of 830-360/min. , and dropped to 810/min. as speed and am- 
plitude increased. The increase of flutter amplitude was 
probably favored by the method of mounting of the aileron 
connecting strut, which sloped about 25° toward the plane 
of the struts. Thus tending of "the upper wing caused mo- 
tion of the lower aileron. 

The flight records ranged "between altitudes of from 
4,000 to 600 m, so as to establish the effect of air densi- 
ty. Flutter- started at 260 to 275 km/h. The reduced fre- 
quency of the lower wing was found at 

' * 0 • 2 

0) = 0.87 (^A (5) 

The air-density effect p is therefore relatively small. 

f )- Cant ile ver biplane gLlA-^ScJiTOlbe^^ This type, 
built since 1927, had been in service quite awhile when, in 
the spring of 1932, several of them developed flutter be- 
low the level top speed tfhich could not be called mild 
because it resulted in fracture of the ailerons. An air- 
plane of this type was therefore subjected to an dscilla- 
t ion t est . , 

The dangerous mode lies at 675/min. It is the symmet- 
rical fundamental bending mode of both wings (fig. 7), The 
nodal line lies far forward of the leading edge of the wing. 
The ailerons are ■ in ' phase opposition; their chord is 240 mm, 
their e.g. is 103 ram behind the hinge line. This results 
in a strong mass coupling between wing bending and aileron 
motion. The mean chord of the extremely oscillating wing 
tips is 1.3 m. The experiments were temporarily inter- 
rupted to permit the airplane to take part in an air cir- 
cus. During this air circus in July 1932, it was stunted 
at speeds up to 200 km/h without developing flutter;- but as 
soon as the pilot started to land, it suddenly began to 
flutter very severely at 145 km/h, which ended in the 
breaking of the ailerons and damage to the plywood covering. 
The flutter continued up to 100 km/h speed. This case 
shows very clearly the unpredictability of flutter. 

The reduced frequency at start of flutter is 



,0).= 0.1885 ~~— "= x ' 14 
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while at the end it reaches the high value of (D e = 1.55, 

g_) Cant ilever_Do _1G jnonq^lane This is an all-metal 

high-wing design. The three- spar wing is traced by short 
struts. The cantilever' length is 72 percent of the semi- 
span. The leading edge of the wing is approximately a 
semiellipse, the trailing edge is straight. During a 
flight on September 9, 1932, it developed such a severe 
case of vicious flutter at about 450 kia/h, that both wing 
tips broke off to the length of one chord. While level- 
ing off from a dive at 2,500 m to 1,500 m, the pilot noted 
oscillations on the ailerons as far as the wing tips which, 
within about 3 seconds, resulted in broken wing tips and 
ailerons. The pilot was able to land safely and the air- 
plane was subsequently repaired and tested on the oscil- 
lating stand. 

The wing reveals a series of oscillation modes in the 
500 to 1,250/nin. frequency zone, whereby the nodal line 
gradually shifts from the front toward the rear spar. Al- 
though the aileron, with a chord of 355 mm, has its e.g. 41 
mm behind the hinge line, these modes do not induce flut- 
ter because the aileron control is very rigid (push rods), 
so that the aileron motion does not build up to large am- 
plitudes at these frequencies, The dangerous mode lies at 
1,400 to 1,500/min. At 1 , 400/min. the nodal line is exact- 
ly coincident with the principal line of failure of the 
wing tips, which slopes 30° outward and backward from the 
leading edge in the direction of flight. At 1,500/min. 
the outer nodal line, in form of a quarter circle about 
the wing tip, is in part coincident with the line of the 
secondary failure. The inner nodal lino runs from the 
point of intersection of the trailing edge and plane of 
struts at an angle of 30° outward, and passes directly 
through a region in which the internal bracing was broken 
(fig. 8), 

Without the lines of failure as clues, it would in- 
deed be difficult in this case to ascertain the dangerous 
mode from the static oscillation test alone. The next 
section contains various factors which should help to fa- 
cilitate this decision. 

Another source of error lies in the estimate of the 
mean chord of that part of the wing which oscillates most 
severely, because of its pronounced taper in plan. Ap- 
praising the mean chord of the severed wing tip at 1.6 m 
and the flutter frequency at 1 , 450/nin . , the reduced fre- 
quency is 
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oo = 0.1885 £i§.Q_X_i._6, = Q - 97 

450 

hlJBr&ce&J^ This airplane is the braced 

high-wing type. The wing structure is of duralumin; has 
two spars, and is covered with fabric. The ailerons are 
split. The aileron chord from hinge line to trailing edge 
is 300 mm; the aileron e.g. is 35 mm behind the hinge line. 
With this e.g. position, vicious flutter started at 290 
kin/h, and persisted to 120 km/h. The rivets of the torsion 
structure were sheared' off from the* outboard aileron hinges 
of both wings. This explains perhaps the rather extended 
range of speeds during which, flutter persists. With per- 
fect mass balance, a speed of 340 km/h had previously been 
obtained without* flutter . 

The airplane was then tested on the oscillating bench. 
Asymmetrical and' symmetrical fundamental bending modes oc- 
curred at 500/min. and 580/min. The natural frequency of 
the ailerons lies at ' 7"30/min, -■ The dangerous mode is the 
torsional oscillation of the wing at 835/min. The nodal 
line runs over the entire span between front and rear spar. 
The aileron amplitudes are high (fig. 9). 

. . For a 1.56 m wing chord, the reduced frequency at in- 
cipient flutter amounts to 

CO = 0.1885 1^56 _ 0<85 

as against the abnormally high u e = 2.05 at its termina- 
tion; the failure of the torsion structure itself may per- 
haps have lowered the flutter frequency. 

il_4CJL2 JE_cant conoj^lane^- This is a cantilever 
high-wing design of wood with tapered wings. While compet- 
ing in the 1932 International Challenge Contest, it devel- 
oped a mild case of flutter "at 220 km/h, but only in rough, 
gusty weather. In fair weather it reached a speed of 270 
km/h without flutter. 

The mean chord of the oscillating wing tip is 1.4 m; 
the aileron chord is 300 mm; the aileron e.g. is 112 mm 
behind the hinge line. 

_ Tested on the oscillation bench, this airplane revealed 
bending oscillations with very indefinite resonance' condi- 
tions at 585 to 830/min. frequencies. The dangerous mode 
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apparently lies between 800 and 830 /rain. , "because then the 
aileron notion caused "by mass coupling and elasticity gf 
the control cables, has a phase difference of about 90 
against the bending oscillation; that is, is in resonance 
with the bending oscillation (fig. 10). T&e reduced fre- 
quency is 

0) = 0.1885 81 5 £ Ui = 0.98. 
220 

ilJDo 12 cant 11 ever.__moncj3lane JjLibel le.!l.- This is a 
high-wing all-metal amphibian. The two- spar wing is ta- 
pered in plan form. 

The airplane showed vicious flutter on September 27, 
1933, at 180 kra/h. The oscillations started when the pi- 
lot opened the throttle after leveling off from a glide. 
The oscillations were so severe that one aileron jumped 
out of its hinges and both wings were badly damaged. The 
wing flutter was preceded by tail buffeting, initiated ap- 
parently when opening the throttle, and which in turn- 
started the wing flutter. The pilot made a safe landing, 
however, after which the airplane was repaired and tested 
on the oscillation bench. 

The dangerous oscillation mode of the wing lies at 
580/min., which at the same time is the principal resonance 
mode of the horizontal tail surfaces. It is an antisym- 
metrical bending oscillation; the nodal line starts at the 
inner aileron and runs outwardly at an angle of 15 in the 
direction of flight (fig. 11). The aileron oscillates in 
torsion. The mean chord of the outer oscillating part of 
the wing was estimated at 1.3 m; the aileron chord is^OO 
mm. The aileron is not mass-balanced; its e.g. position 
was estimated at 100 mm behind the hinge line. The aile- 
ron control cables are not very rigid. The reduced fre- 
quency is 

0) = 0.1885 m.JL±t* = 0.79.- 

^ 1 _M _2 8 _m.o nojpl a ne . - This is a cantilever low-wing de- 
sign', of "duralumin with wings tapering in plan only. 

After extensive testing, the airplane developed a mild 
case of flutter at 220 km/h, which was started by the 
(mass-) unbalanced ailerons. It stopped when the speed 
was reduced to 180 km/h. The pilot had the impression 
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that the flutter started each time after a "bump, even, when 
very slight. The engine r.p.in. was 1,750 at the "beginning. 
Later the airplane was dived to 250 km/h at 1,950 r.p.m.' 
without developing flutter. 

The. airplane was then subjected to an oscillation 
test. The wings have a symmetrical fundamental flexural 
oscillation at 480/min. and ant i symmet ri cal "bending oscil- 
lations at 850 and 770/min. The resonance points are very 
clearly expressed. The dangerous mode ^eems to lie at 
770/min., "because this is the frequency at which out-of- ■ 
phase oscillation of the ailerons is first noticed, which 
likens it much to the dangerous mode of the Do 12." The 
nodal line runs outward from a point near the inboard end 
of the aileron at an angle of 10° in the direction of flight 
(fig. 12). 

The mean chord of the outer oscillating part of the 
wing is taken at 1.05 m, the aileron chord at 350 mm; the 
e.g. is 140 mm behind the hinge line. The reduced fre- 
quency is 



It is planned to make flutter measurements in flight on 
this airplane in order to determine the flutter frequen- 
cies exactly. 

Ll^iE.lane_S ,.2 4_^Kiebit z^.- This is a biplane of wood 
construction, braced in one plane. In the spring of 1932, 
the airplane went into a long, unexpected dive with' a burn- 
ing engine and started to flutter, finally breaking the 
cellule. The calculated terminal velocity is 280 km/h. 

Another airplane of the same type was subjected to an 
oscillation test. At 490 and 515/min. the whole cellule 
started to oscillate; at 825 and 1000/min., the overhang 
went into flexural oscillations. The dangerous mode lies 
at 1215, because this was the frequency at which the aile- 
rons first revealed phase opposition. The nodal line runs 
from the intersection of the strut plane and trailing edge 
to the first third of the edge strip (fig. 13). 

The mean wing chord is 1.18 m, the aileron chord 240 
mm, the e.g. of the aileron system is 21 mm behind the 
hinge line. The reduced frequency is 



oo = :0.1885 



770 X_1^05 
220 



0.69. 



= 0.1885 



!215_X_1.18 

• 280 



= 0.97. 
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Biplane Ar . 6 6 C.- This is a traced "biplane with a 
smaller lower wing. In a dive at 340 km/h,,the airplane 
started to flutter, which. led" to the beginning' of a frac- 
ture' oi : the lower front spar as well as of the plywood 
covering on the lower side of the wing. Since the stick 
did not oscillate, the mode was symmetrical. The chord of 
the lower wing is 1,65 m;. the e.g. of the' aileron lies 45 
mm behind the hinge line. 

The airplane was subjected' to an oscillation test by 
the manufacturer. The bending oscillation of the lower 
wing at 790/min. was considered as the dangerous- mode. At 
this mode, the inboard part of the wing pivots roughly 
about the. front spar, while the overhang bends. 

The reduced frequency is , 
(0 = 0.1885 mJL-UJ®. = 0.73 

After finishing this report, flutter was again observed 
after the ailerons had been completely mass- balanced and 
were perfectly quiet in the oscillation test. Following 
some minor changes in the shape of the' aileron, it devel- 
oped vicious flutter at 420 km/h, leading to complete de- 
struction of the cellule. This gives 

(i) = 0.1885 £££ JL U15 = 0#58 

This might have been a' case of flutter in combined bending 
and torsion, although not without some probable aerodynamic 
coupling effect of the aileron motion. 
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TABLE II, Observed Cases 
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TABLE Ila. Observed Cases 





lype 


Mode of 
lat ion 


Rigidity 

01 

cont rols 


A spect 


Wing 
struc- 
ture 


1 


DP 9 


symmet- 
rical 


small 


vicious 


wood 


2 


He 8 


ir 


II 


n 


n 


3 


L 78 


ant i sym~ 
met rical 


great 


mild 


1! 


4 


He 60 


- 


- 


vicious 


metal 


5 


He 46 c 


symmet- 
rical 


small 


mild 


wood 


6 


KL 1 A 


H 


ii 


vie ions 


ii 


7 


Tin 1 0 


an t x s j hi 
met rical 


g r eat 


ii 


met al 


8 


L 102 


- 


small 


it 


ii 


9 


AC 12 E 


symmet- 
rical 


M 


mild 


wood 


10 


Do 12 


ant i sym- 
metrical 


tl 


vicious 


metal 


11 


U 28 


it 


great 


mild 


ii 


12 


S 24 




small 


vicious 


wood 


13 


Ar 66 C 


symmet- 
rical 


1! 


ii 


n 



2. Conclusions 

The numerical data of these 13 test cases are append- 
ed in tahles II and Ila* The mean value for incipient 
flutter is, according to table II: 

03 m = 0.90 ±0.12, 



1\ T .A.C.A. Technical Memorandum No. 782 



25 



The majority of cases are vicious; only 30 percent 
are mild (table Ila). The structural material does n6t 
seem to have any effect. 

When the stiffness of the ailerons and their controls 
is great, flutter occurs in an anti symmetrical mode he- 
cause only then can the ailerons oscillate freely and 
transmit energy. 

If the stiffness of the aileron' control is small, 
flutter may also occur in a symmetrical mode, "but then on- 
ly at a frequency high enough above the natural aileron 
frequency to permit motion. in phase opposition. Oscilla^ 
tions of the whole cellule of a biplane about the vertical 
and the transverse axis may not necessarily .lead to flutter, 
even when combined with additive wing torsion. The danger- 
ous mode, however, is frequently the first natural mode 
of the wing independent of the cellule. 

The reason for this behavior lies. in the damping of- 
the oscillations through the precessional moment of the 
rotating propeller. Heretofore, all airplane o scillat ion 
tests have been, almost without exception, static tests, 
i.e., with the engine standing still. Under these condi- 
tions, a number of cellule * oscillation modes may develop, 
during which the fuselage oscillates slightly in torsion. 
In flight with full r. p.m., the precessional moment of the 
propeller damps such oscillations very effectively, so 
that flutter is very rare. The KL 1 A revealed such cel- 
lule oscillations on the stand at the dangerous frequency. 
The sudden entry of flutter when starting to land, was 
probably attributable to the diminished damping of the pro- 
peller as a result of the smaller r.p.m. One two-engine 
airplane tested on the oscillation bench manifested a 
marked difference in oscillation modes, depending on wheth- 
er the engines were running or not. On the other hand, 
purely symmetrical wing oscillations, during which the 
thrust axis is merely translatory, do not prevent flutter 
as proved by the Ar 66 C airplane. 

Table II also gives the natural frequencies of the ai- 
. leron^ oscillations n for various airplanes , together with 
the width of the resonance curve An defining the damping 
of the flutter oscillation. This width is measured at 71 
percent of the maximum amplitude. 
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IV. DETERMINATION OF LIMITING- CONDITIONS FOE FLUTTER 



Section II shows that further analytical. treatment of 
concrete cases of flutter is impractical "because of the 
many secondary circumstances which are not amenable to nu- 
merical treatment. Even so, it is useful from the point 
of view of flutter prevention, to establish analytically 
the limiting conditions under which flutter. is possible. 

. Flutter is obviously possible only when the oscillat- 
ing- -wing ' is able to take energy out of the air stream in 
order to . equalize the ever-present damping losses. For a 
general survey, the calculation may be restricted to the 
limiting conditions of energy absorption in two-dimensional 
flow. 

The reduction of lift due to tip vortices of the os- 
cillating wing tip, will have the effect of increasing 
the damping and narrowing the range in which flutter is 
possible. This additional damping effect can be estimated, 
although it is neglected in the following derivation. 

' - » * 

* 1. Notation 

£ I m stroke ( tending) amplitude of (aerodynanic) 

neutral axis ( quart er-chord, point) . 

B/ . torsional amplitude of wing. 

C deflection amplitude of aileron. 

Id m effective span of oscillating wing tip, 

aileron^ span, 

damping factor of wing. 

damping factor of aileron. 

torsional stiffness of wing at abscissa y. 

distance from nod^l line of wing to center 
line of wing. 

energy of oscillating airplane, 
energy of wing. 



m 



D (y) 



kg m 2 
in 



E 



m kg 
m kg 





"KT ft 

Jn « A. 
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m kg 


energy of aileron. 


P 


m s 


wing area. 


g 




damping phase angle of wing oscillation. 


h 




damping phase angle of aileron oscillation. 


k 


1) 

2) m 


stiffness constant. 

radius of gyration of wing section, referred 
to nodal line. 


k E 


m 


radius of gyration of aileron, referred to 
aileron nodal line. 


I 


m 


semichord of wing. 


l x .. 


. • I g 


energy co efficient s . 




m kg/ s 


mean aerodynamic energy. 


L dm 


n kg/ s 


mean damping power. 


M d 


in kg 


torsional moment. 


m 


kg s s /m 


mass, 


m F 


kg s p /m 


mass of oscillating wing tip. 


m R 


kg s s /m 


aileron mass. 


n 


1/min. 


oscillation frequency. 


n o 


l / mm . 


oscillation frequency at resonance. 


An 


l/ min. 


width of resonance curve in oscillation 
test. 


P 


m 


absolute oscillation amplitude. 


s 




position of wing e.g. "behind neutral point 
(quarter- chord point^. 


S R 




position of aileron e.g. hehind hinge line. 


t 


ra 


wing chord. 




m 


aileron chord. 
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v m/s; km/h flight speed. 

y m coordinate along span* 

£ index of amplitude ratio for wing. 

£ R index of amplitude ratio for aileron. 

\x = ~- leg/m 3 specific weight of wing. 
Ft rn 

v l/s frequency of wing oscillation in radians. 

p leg s p /m 4 air density. 

p o air density, at sea level. 

T = ^ aileron chord ratio. 

cp l) angle of twist. 

2) phase angle between tending and torsional 
o sc illat ion. 

^ phase angle "between "bending and torsional 

oscillation of aileron. 

X phase angle "between torsional oscillation of 

wing and torsional oscillation of aileron. 

$12 functions of aileron chord ratio. 

0) reduced' frequency. 

00^ maximum value of reduced frequency. 



2. Aerodynamic Energy of Oscillation 

With A I denoting the stroke (bending) amplitude, 
B and C, the amplitudes of angular motion of wing and aile- 
ron; and index 1 signifying the real part, index " the im- 
aginary part of these amplitudes; a "bar - their absolute 
magnitude, the^tirae average of the aerodynamic energy of ■ 
oscillation is, according to (16) and (17): 



*The theoretical development of the ensuing formulas is to 
he published in a separate report. 
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L m = t 71 P v v * 1 10 & ll + B 3 + ° ts 

+ (A«B« + A" B" ) 1 4 + (A'B" - A"B«) l s 

+ (A'C + A" C") l e + (A'C" - A" C 1 ) l 7 

+ (B'C + B" C" ) l 8 + (B'C" - B" C" ) l 9 

where the energy coefficients: 

li = 1 + T" 

l 3 = 1 

I 



= _1 
2rr 



CM. £ + 



(i + T t) + *xJ 



U s= + 2 + T I 



, 5 = l_t_Tl _ T „ 

5 0) 



l s = 1 [$ x SI + ^_+le (i + t«) + * 3 ] 

T»] 



l 7 = 1 C * x 14-1 



2 



(6) 



(7) 



Tho functions T ! (o)). T"(co) "based on Hankel ! s cylin- 
der functions, are given in tahle III. The functions 
$n( T ) derived from trigonometrical functions, are given 
in tahle IV. Prom these the energy coefficients l n have 
"been computed for aileron chord ratios of T = 0.15, 0.20, 
0.25 with a 20-inch slide rule (tahle V) . The energy is 
negative when taken out of the air stream. 
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TABLE III. Functions T'(o)) » T "(oo) 



00 


0 


0 .2 


0 4 


0 fi 


\J . o 


3 n 
1 t u 




T ' 


1 


0 .455 


0 .250 


0.158 


0.108 


0.079 




J II 


0 


- .377 


- .330 


-.276 


-.233 


- .201 




W 
T 


1 


1.2 

0 .0 60 


1.4 
0.047 


1.6 
0.038 


1.8 
0.031 


2.0 

0.026 


oo 

0 


rp Tl 


-.175 


- .156 


-.140 


-.126 


-.115 

, — 


0 



TABLE IV. Junction § n (T) 



T 


0.10 


0.15 


0.20 


0.25 


0 .30 


*i 


1 .244 


1 . 510 


1.727 


1.913 


2,076 




.334 


.610 


.935 


1.299 


1. 698 




.164 


.296 


.447 


.614 


.793 


& 4 


.0264 


.0719 


.1459 


.2518 


.3924 


$S 


1.080 


1.214 


1.280 


1.299 


1.283 


$6 


.615 


1.077 


1.577 


2.094 


2.612 


$ 7 


.0506 


.1347 


.2672 


.4507 


. 68 60 


$0 


.0079 


.0192 


.0400 


.0707 


.1129 


$9 


.0459 


.1246 


.2519 


.4330 


.6718 




.0546 


.1803 


.4180 


.7978 


1.345 




.0056 


.0281 


.0876 


.2128 


.4371 



N.A.C.A. Technical Memorandum lib. 782 31 



TABLE V. Energy Coefficients 



00 




0-4 


0 . fi 


0 9, 


1 0 
x . u 


1 P 

X . 


J. , D 






7 

0 1 


1.2 50 


1 .158 


1 . 108 


1.079 


1 .0 60 


1.038 


1 .026 






1 


X 




X 


X 


X 


X 






1 #425 


1 . 698 


1 . 817 


1.878 


1 . 914 


1 . 951 


1.968 


— 


^5 


3. 455 


2 .205 


1 . 618 


1.279 


1 .059 


. 788 


. 628 




u 3 


0 .00829 


0 .00880 


0 .00904 


0.00916 


0 .00923 


0 .00931 


0 .00935 


LO 


l« 

b 


-•1773 


- .010 6 


.0652 


.1059 


.1300 


.1561 


.1696 


rH 
* 

o 


l 7 


1. 534 


.954 


. 688 . 


. 538 


. 442 


. 326 


.258 




o 


.1926 


. 1934 


.1938 


. 1939 


. 1941 


.1942 


.1943 




*9 


. 37 57 


.3796. 


. 3815 


. 3825 


. 3832 


. 38 40 


. 3846 




" 3 


0.0195 


0 .0207 . 


0 .0212 


0 .0215 


0 .0217 


0 .0219 


0 .0219 




1 „ 


- .1175 


.0 695 


.1542 


.1993 


-2265 

. (*J £j U \-J 


.2553 


2697 


o 

0) 

d 


l„ 

7 


1 .769 


1 . 104 


, 798 


.62 4 


. j. 0 


.379 


- 300 






.2 940 


.2956 


.2964 


.29 67 


-2970 


.2973 


-2976 




l 9 


.3854 


.3934 


. 3970 


.3993 


.400 6 


.4024 


.4034 




^3 


0.0377 


0.0400 


0 .0410 


0.0416 


0 .0419 


0 .0422 


0 .0424 


0.25 


Is 


-.0343 


.1682 


.2596 


.3078 


.3374 


.3686 


.3841 




1.975 


1.235 


.895 


.701 


.576 


.42 6 


.335 




I's 


.407 


.410 


. 412 


.412 


.413 


.413 


. 41 4 




*9 


.3746 ' 


.3886 


.3952 


,3991 


.4016 


.40 45 


.40 6 
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-3. Damping 

As the damping energy was not measured directly in 
the oscillation tests, the energy absorbed "by material; 
damping and external friction must he estimated from the 
phase difference "between elast ic - f or ce "and deformation; 

The previous assumption of a constant phase differ- 
ence, gives an elliptic hysteresis loop independent of fre- 
quency; which Agrees closely with observed facts. The' 
product of this phase difference and tt equals the usual 
logarithmic decrement. 1 A detailed exposition of these 're- 
lations will.be found in. .B. v. Schlippe f s article (refer- 
ence 25) . 

• The damping energy is obtained from the relation 

tt g ™nean~energy of oscillating system 

• ** * 

The phase angles g and h of the wing and pf the aileron 
oscillations may be determined "from the width of the reso- 
nance curves of the oscillation test. If the width of the 
resonance curve An is measured at 71 percent of the reso- 
nance amplitude (fig. 14), the phase angle is 

g = A ^ . ■ (-8) 

At higher amplitudes tlie phase'angle increases, the in- 
crease being, as a rule, greater for wooden wings than ^ 
for met'al wings. 

An exception to this is the damping due to friction 
in the aileron hinge .bear ings-. With constant fri'ctional 
moment 'the damping is proportional to the (angular) ampli- 
tude C, while the energy of -the oscillating aileron rises 
proportionally "to C 2 . 

By the same argument, the 'phase angle h of the ai- 
leron .is "inversely proportional to -the amplitude. Thus, 
at very low amplitudes, the phase angle of the damping can 
be so great as to make flutter impossible, Flutter cannot 
set in unless some outside cause imparts a momentum of a 
certain minimum rise to the aileron. At high amplitudes 
only the residual damping due to material damping of the 
control cables remains. With careful installation of the 
ailerons and cables this effect is less pronounced. 
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Exactly corresponding averages must tie formed for the am- 
plitudes -B and C. The mean damping energy is 

L dm = v{s V h E ) = ^ tt p V 2 v I 3 (0(d B S +d C 2 ) (12) 



For the sake of "brevity, the following nondimensional damp- 
ing factors have been introduced into equation (12): 



d p = g 



[tt p I* 



dp = n ] 

K TT p I 



. 1. . 
— + — + — s 

h 8 I 



* 4 tt 2 "b 



(13) 



Table VI gives damping factors for several airplanes. 
The factor dp 1 applies to pure "bending oscillation. The 

torsional oscillation about the neut fal po int is estimated 
at 0»4 dji 1 "because the nass distribution along the chord 

is not very unlike in the various types. The. aero dynamic- 
ally effective wing span b may bo estimated from the 
dangerous oscillation mode. Only a portion of the outer 
oscillating part of the wing "between wing tip and nodal 
line or plane of struts can he considered as aerodynanically 
effective, "because the trailing vortices lower the efficien- 
cy of the wing tips as compared with two-dimensional flow. 



It should he emphasized that this estimate of the 
damping is quite rough and therefore merely affords the ap- 
proximate magnitude of the damping energy. For this rea- 
son, the direct determination of the energy of damping 
from oscillation tests is very much desired. 

4, Ratio Qf Amplitudes 

In order to he able to compare the amplitude ratios 
of the forced oscillations with those of the oscillation 
test, the elastic o sc illation of a flat plate covered with 
a mass was investigated. The plate is assumed to.be piv- 
oted about the axis P and elastically restrained against 
torsion, so that it oscillates with a natural torsional 
frequency n (fig. 15). 
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The phase angles of wing "bending and wing torsion 
damping may he considered as being about equal. Table VI 
gives the values of phase angle g obtained from oscilla- 
tion tests with constant excitation and with impulse ex- 
citation (deflection tests). They are multiples of the 
values measured on smooth test bars, because in built-up 
members develop additional losses at the joints* On the 
average g = 0.061. The phase angle of the aileron damp- 
ing is according to (6) h = 0.20, and according to a res- 
onance curve, h = 0.08. 

The mean energy content of the oscillating airplane 
at the dangerous mode may be determined by forming the in- 
tegral 



E 



2 J 



dm 



(9) 



over" the whole airplane, whereby 
amplitude of the mass element dm 
lation test. This calculation is 
that in most cases an estimate is 



p denotes the absolute 
measured in the oscil- 
so tedious, however, 
preferable. 



Visualizing the outer end of the wing as a rigid 
plate oscillating in torsion at amplitude B about its 
nodal line while the aileron oscillates in torsion about 
the hinge line at amplitude C, the energy content of wing 
and aileron is: 



% = f- U k 3 '+ TT P V 



E H = ~ [m R k R 2 + p I 4 



* / 



(10) 



If the wing oscillates in pure bending, then k = e = °° ? 
B k = B e = A. 

Very often the oscillation amplitude and the wing 

chord are very variable near the tip. As the energy trans 

fer involves the square of the amplitude, it is advisable 
in this case to form the mean values: 
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The three following limiting cases oust "be considered: 



1« n 0 = 0 , free pivoting about the axis P. 
~ = phase angle 9 = rr. 



B 

2 , n 0 = n , 
A _ 



forced frequency = natural torsion. 
" frequency. 



- = g, phase angle cp= ^ . 



B 



♦very high, great torsional stiffness. 
& ™- J 1 + g s » phase angle cp - g. 



n c 



Similar relations hold for a plate with attached aile- 
ron, which is pivoted with elastic restraint about the 
hinge axis. With A R I (fig. 16) as the forced amplitude 
of the aileron hinge, the amplitude ratio "of the forced 
oscillation is: 



C 



(15) 



The values for the amplitude ratio in the three limit- 
ing cases correspond to those given above. The cases be- 
tween the first and second limiting cases are of particular 
significance for the forced oscillation, because here the 
phase angle lies in the second quadrant; that is, it ap- 
proaches the optimum phase angle which almost always lies 
in the third quadrant. This is borne out by experience. 
Flutter usually occurs either antisymmetrically - that is, 
with freely oscillating ailerons, so to say - or symmet- 
rically, at a frequency which lies above the natural oscil- 
lation frequency of the ailerons* 



The values of the amplitude ratio obtained from the 
oscillation test are, on the whole, smaller than according 
to limiting case 1, because the elasticity of the control 
system shifts the oscillation mode toward limiting case 2. 
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No. 


Type 


g 


kg/n 3 


/ Oi t_Li \ 




— 




1 


, — . 

DP 9 


0 »0 oU 


« Q 


U • DO . 






2 


He 8 


• 030 


<3 • 4 


T 


1 A 




3 


L 78 


.034 


6 . y 




T R 
• ID 

• 




4 


He 60 


ACT 

• 0 51 


c5 • b 


AO 




\J 9 \J \J £j 17 


5 


He 46 


.047 


f ■ <o 


• r c, 






6 


ZL X A 


• 0 67 . 


b # o 


OA 


• x t> 




I 


D 0 10 


A AT 


■ftp 


• op 


.♦19 
* x ^ 


.0017 


8 


L 102 


.042 


5.0 


.46 


.19 




9 


AC 12 E 


.147 


5.4 


1 i 67 


.21 


.0046 


10 


Do 12 


.086 


5.1 


.96 


.31 




11 


M 28 


.078 


5.4 


.92 


.33 


.0134 


12 


S 24 


.0 58 


4.0 


.52 


.20 


.00 62 


13 


Ar 66 C 













Assume that axis P executes forced vertical oscilla- 
tions at amplitude A I and frequency n. , The, amplitude 
ratio of this forced oscillation is: 



4 = 

B 



. G-i- ^ - *) 



& = — 



n k 3 , + |.TT p I * 



-, 1 i* v 

Q S I + 5" TT p I o 



(14) 



The comparative factor $ ' depends, only on the mass 
distribution and on the air density. The second terms 
represent the effect of the co-oscillating air masses* 
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It should be noted that finite values of 



R 



are 



still obtainable even if the aileron c.g* lies in the 
..hinge line, "because the co-oscillating air nasses exert a 
.mass— coupling moment. This explains the fact that an ai- 
leron, whose static e.g. lies only very little "behind the 
hinge line- can still excite the. o sc illations to flutter 
(cf. section II 1 , 1 , g,h, c) . 

For the customary method of estimating the critical 
velocities, as well as for the new method using oscilla- 
tion tests and recommended hereinafter, it is important to 
know, whether flutter occurs in approximately the same mode 
as the oscillation of equal frequency on the stand. 

A very simple way is to visualize the wing tip as a 
mass-elastic system. Posing the torsion of the wing tip 
about the nodal line at 



P -Be 



ivt 



.whereby amplitude B is real, the inertia force moment is 



dt 



^ = ~ 9 B v 



s ivt 
e 



the elastic moment 



f 0 e ig Mt 0 B e iUt+iS , 



and the ezcitational moment 



= M e 



iut 



In steady oscillation, the sum of these three moments 
must always disappear; that is, 

- m B + f 0 B e le + U = 0 ' 

Separating this equation into real and imaginary parts, 
affords 

- e b V s ; + f 0 b cos g + in = o 

i f 0 B sin g + i M" = 0 

How, the aerodynamic exciattion could occur in such 
a manner that the "blind component" !J» = 0. Then the 
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oscillation frequency with aerodynamic excitation would not 
differ at all from the natural frequency and fron the fre- 
quency H ! = 0 of the oscillation test. In actual flutter, 
it usually is H 1 < 0. The phase angle most favorable for 
the excitation lies Generally in. the third quadrant, as 
will he shown below. The oscillation frequency is, as a 
result, lower. On the other hand, the frequency increases 
with the flying height because the air nass co-oscillating 
with the wing ,. decreases with decreasing density. 

The estimate will he fairly correct when assuming ^ the 
phase angle of the exciting moment at cp ~ 225 ; that is, 
posing 

Hi == M" — f 0 A sin g r 

the average for 12 wings is sin g ~ 0,06. The ^ change 
of frequency in this case is -3 percent; that is, still 
within the resonance width An. Model experiments and 
flight measurements in several instances proved that ^ these 
assumptions agree with the facts in approximate magnitude, 
particularly for coupled wing and aileron oscillations* 

Greater blind components and consequently greater dis- 
crepancies in frequency are to he expected when the damp- 
ing is great and the phase angle of excitation is close to 
180° or 360°. This is the case, e.g., on model wings with- 
out ailerons which are not dynamically similar; have high 
frictional damping and oscillate with a low reduced fre- 
quency. The limiting case is the "oscillation" with 0 
frequency - the aperiodic twisting off of the wing due to 
static torsional instability. In this case, 

Hi = - f Q B cos g. 

However, the wing flutter, observed up to date, has 
been bo far from this limiting case that the assumption of 
identical modes for actual flutter and for oscillation 
test may be considered as, a close approach, provided, of 
course, the oscillation test is made with different meth- 
ods of excitation since, for example, symmetrical modes 
appear only indistinctly or not at all, with ant i symme t r i- 
cal excitation. Furthermore, it should not be expected 
that oscillation modes due to aerodynamic coupling, are 
reproducible on the o s c illat ion bench. 
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5. Criterion for Flutter 



Wing flutter with constant amplitude occurs when the 
energy, of damping equals the input of aerodynamic energy; 
that is, when 



L m + L dm = 0 



(16) 



If, for simplification, we set A ! = A, A ,r - 0, we ohtain 
from equations (6), (12), (16) the condition under which 
oscillations with constant amplitude occur; 

A S (l x + oo d F <) 4- B S (l 2 + 0) 0.4 d F M. + 0 2 (l 3 + 00 d R ) 
+ A B 1 I4 + I B" I s + A C 1 l G + A C" 1^ 

+ ( B 1 C + B» 0") l„ + (B< C" ~ B» C) l 9 = 0 (17) 

As the oscillation test is primarily a means for re- 
cording the oscillation amplitudes, it is advisahle to 
consider the amplitudes B and C as independent variables 
and to set 



B 1 = B cos cp 
B" - B sin cp 



C T = C cos \|/ 
0" = C sin t 



The phase angles cp and \|/ should "be so determined 
that the energy assumes an extreme Value,. Partial differ- 
entiation of (10) according to cp and f gives: 



(18) 



A(~l 4 sin cp 0 + l s cos q> 0 >+ 0 

[l a sinty 0 ~9 0 )- l 9 coa(t 0 -q> 0 )] = 0 
A(- l G sin \|/ 0 +1 7 cos t 0 ) " B 

[1, sin(y 0 -cp 0 )- l 9 cos(t 0 -<P 0 )^ = 0 

From (17) and (18) the limiting values of the ampli- 
tudes B and 0 may he computed as functions of the reduced 
frequency. Admittedly, the calculation is quite complicat- 
ed for three degrees of freedom. 

Restricted to two degrees of freedom, the results are 
more simple and elucidating* This limitation is particu- 
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larly permissible when the aileron is the main cause of 
flutter and the rest of the wing simply oscillates with it, 
Then the ensuing mode of wing oscillation may be said to 
have only one degree of freedom. The wing then oscillates 
about some fixed nodal line as in the oscillation test. 
This has also been observed in model tests. The four fol- 
lowing modes of oscillation with two degrees of freedom, 
are investigated. 



1. 


C = 


0, 


wing bending + wing torsion. 


2. 


33 = 


0, 


wing "bending + aileron motion. 


3. 


A = 


0, 


wing torsion about quart er- chord point 
+ aileron motion. 


4. 


I = 




B ' , B" = 0, wing torsion about three- 



quarter chord point + aile- 
ron motion. 

Writing the extreme conditions (18) in (17) results 



in: 



1. A 3 (l a + 00 dp' )+B \l s + 0) 0.4 d F » )-A-bJi^ +l 5 2 = 0 



sin cp 0 = - —=====., cos <p 0 



2. V{\ + ood F ')+C 2 ('I Q + 0)d R )-AC yig +lf = 0 

I l- 
sin % = - -=r-JL~—-, cos ^ 0 = 7—7— 



3. B (I a +W0.4 d F ') + C (I3+00 d R )-BC A /l 8 3 + I 9 s = 0 

4. B 3 ( I W 1.4 dp') + C 2 (1 3 + to d R ) 



- BC/Og - l 6 ) + (l 9 
sin X 0 - 



l 7 ) = 0 



*8 - l e 



/u a - i G f + (i B - i 7 ) 2 



cos 



a * e' 

I g " I . 



> (19) 
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Excepting l 6 -,th.e energy coefficients are consist- 
ently positive; therefore, the optimum phasp angles usu- 
ally lie in. the third quadrant. Tor mode 2, the optimum 
phase angle passes from the second into the third quad- 
rant at a) ~ 0*5; , for the fourth mode, it passes from the 
fourth into the third quadrant at 03= 1.5. 

6 . Energy Graphs 

In order to circumscribe the influence of damping on 
flutter, the calculation is "based on the following extreme 
value s : 

dp= 0 0.25 1.25 

*E = 0 .0025 .025 

The values d = 0 are only of theoretical interest. 
The corresponding curves show the range within which 
energy may he taken from the air stream. 

The four oscillation modes finally afford energy 
graphs, with four curves each, for the different damping 
values. The ahscissa'is the "logarithm of the amplitude 
ratio; the ordinate is the reciprocal value of the reduced 
frequency , 

1 = JL 

00 vl 

In the following, the energy graphs (figs. 17-23) are 
compared with experience. 

a L 0 s cilia t_i o n. :__w 1 S&_h en d i ng,_an d _w i ng_t o r s i o n . - Fig- 
ure 17 shows that the hest condition for flutter exists at 
the amplitude ratio' a/b = 0.9. This corresponds approxi- 
mately to a pure elastic oscillation ahout the three-quar- 
ter chord point. The reduced frequency a> h = 0.85 is not 
exceeded even with very low damping dp' = 0.35, and may 
therefore he considered as the practical limit. With 
greater damping, as it occurs particularly on models, the 
upper limit drops to 0) = 0.52, so that in especially un- 
favorable cases, flutter. in tending and torsion is likely 
to occur. In the region «" = 0.52 to 0.85, this type of 
flutter has heen analytically investigated on the Junkers " 
A 20 (reference 17) # The tapered wing has a chord of 2.36 
m at the root and a mean chord of 1.70 m at the severely 
oscillating tip. Referred to the wing chord of 1.70 m, 
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the calculation gave in second approach, the reduced fre- 
quency oo = 0.45. The torsional stiffness of the thick 
cantilever Junkers wing" is so great that the critical 
speed would still lie above the operating range even if 
the reduced frequency were twice as high. 

hl_0_s_cillat ion_mod 
This mode" is especially predominant in the He 46, KL U, 
AC 12E, Do 12, M 28. 

Figures 18 and 19 are calculated for aileron chords 
of T = 0.15 and 0,25, and show first that the flutter 
zone becomes smaller as the aileron chord becomes greater. 
According to (15) the backward e.g. position of the aile- 
ron affects the amplitude ratio. The reduced frequency 
is so much higher as the amplitude ratio is smaller and 
the aileron e.g. moves backward, although increases only 
to 0) ~ 2 with very small damping and 0) ~ 1 with very 
great damping. Practically only the right-hand parts of 
the curves come in question, because the amplitude ratio 
of the flutter oscillation cannot drop below a certain 
value determined by the mass distribution of an unbalanced 
aileron. The value A/C = 0.2 to 0.3 may be regarded as 
the practical limit. 

If the amplitude ratio is low the aileron damping ex- 
erts a profound influence on the reduced frequency. At 
A/c = 0.25, e.g., CO increases from 1.13 to 1.4 when the 
aileron damping drops from an initially high value of 
0.025 to the low value of 0.0025. Flutter therefore con- 
tinues until the' speed has dropped to 80 percent of its 
initial value. 

A comparison of the & values as computed for several 
airplanes with the energy graphs, discloses CO values 
which are of the sane order of magnitude as the observed 
ones. As the interpolation for different aileron chords 
and dampings by means of figures 18 and 19 is not very ex- 
act, a closer agreement can be obtained when establishing 
a special energy graph for each airplane, using the damp- 
ing values and amplitude ratios obtained in the oscilla- 
tion test and taking the oo values from these graphs. The 
essential factor is the correct reproduction of the damp- 
ing effect and the backward e.g. position of the aileron 
on- the flutter phenomenon. 
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c) 0 s c i 1 1 a t io n _rno de_:_ W ing _tor s i o n a "bout a guar t er- 

chord point and aileron motion.- This mode occurs appro x- 
imately'~in models' DP "9 , He 87~~He 60, L 102, S 24. The 
fact that 80 percent of these types of airplanes were com- 
pletely destroyed, proves that this node is "by far the 
most dangerous* The energy graphs (figs. 20 and 21) show 
that extremely high values of the reduced frequency are 
obtainable with small damping. If the aileron damping 
drops from an initially high value to a small final value 
after flutter starts, a large excess of energy is availa- 
ble, which inevitably must lead to complete failure. 

One interesting feature is that such flutter is possi- 
ble only within a limited range of the amplitude ratio. 
The practical upper limit is for the aileron chords: 

T = 0*15 B/C 1 0.37 

T = .25 B/C < .45 

This would stipulate a certain minimum distance of back- 
ward e.g. position of the aileron. 

An .aileron with the rather convent ional characteris- 
tics: t = 0.15; k s = 0.12 I ; s R = 0 .0 5 I ; * m R = 0 • 2 5 tt 

p I bq has, for example, according to (15), the ampli- 
tude ratip 



if freely oscillating. According to figure 20, this value 
lies exactly at the point of minimum critical velocity for 
small wing damping, and gives u> = 1.25. In the neighbor- 
hood of amplitude ratio B/ Q = 0.2, reduced frequency 
values up to OJ - 0 • 9 are possible even with maximum 
damping. With small damping, very high co values are pos- 
sible, as actually observed on the L 102 at the end of 
f lutt er . 

d) Osc illation jniojie :_J[in^ 

chord point + aileron motion.- In this r:.o d e no reduced 
frequency values in excess of 0) ~ 1.0 are possible, even 
with small damping (figs. 22 and 23). This mode is there- 
fore less dangerous. 

The aileron chosen above as example, requires twice 
the unbalance ( s R ~ 0.12 I) to give an amplitude ratio 
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B « JL.*B^ .0.89 
C 0.2 C. 

which, in figure . 22, meets the right-hand ; side of the curve 
for small damping at. CO ,=■ 0.62. With the. aileron c.g# 
farther forward , this value could "be obtained only when the 
amplitude C is 'increased at the sane time "by resonance 
of the control, system. 

7. Application to a Practical Example 

tfhese energy graphs are supposed to give a general 
view of %he new method of estimating the critical speed. 
In actual cases , it usually ' involve s oscillations at which 
the nodal lines do not exactly correspond with the third 
nor with the fourth type of the illustrated examples. The 
aileron chords also vary within a great range. However, 
with the aid of the figures given in tables III and IV, 
supplemented "by the damping energy from the oscillation 
test, a particular 'energy graph can be obtained for each 
individual , case . ,. 

Tor example, if the oscillation test shows that the 
nodal line in the outer part of the wing lies at three- 
quarter chord point, that the (angular) amplitude C of 
the aileron motion equals twice the wing torsion amplitude 
B , and that the wing damping factor is dp 1 = 0.25, then 
figure 23 gives for b/c = 0.5 



i- = 1.06 to 1.10 

00 

0) = 0.94 to 0.91 

depending on the value of the aileron damping d R . In this 

case the effect of d^ on the reduced frequency is quite 
small* In fact, the. reduced frequency will be lower than 
the values found from the graph, because the assumption of 
optimum phase angle is not exactly fulfilled. Even so> it 
is possible to estimate the lowest critical speed 'at the - 
observed mode of oscillation, which may prove very valua- 
ble under /certain circumstances . 

Preventative measures against flutter, particularly 
mass balancing of the ailerons, show their effect, in the 
low aileron amplitude C in the oscillation test. the am- 
plitude ratio . B/C can then increase quite easily beyond 
1, so that ho flutter at all is possible for. a nodal line 
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at the quarter-chord point (figs. 20 and 21) or at any 

rate that a substantially higher value of l/oo will he 

found on the right-hand branch of the curves in figures 
18, 19, 22, and 23. 

V. PREVENTION OP FLUTTER 



The previously described conditions for oscillation 
modes of amplitude ratio and damping - while admittedly 
necessary - do not, however, constitute adequate conditions 
for flutter. The thus- est imat ed critical speeds are there- 
fore on the safe side. In the energy calculation , the prem- 
ises were optimum phase angles which, as is known, are 
.closely approached in many cases in practice, especially 
with aileron oscillations. On the other hand, it. is con- 
ceivable that cases may occur wherein the phase angle can- 
not even' approach the optimum value and. in which no flut- 
ter is at all possible, even if the oscillation modes as 
recorded on the oscillation bench, were indicative of flut- 
ter. 

The very simple form of the energy method compared 
with the exact .method , is simply the result of omitting 
the elastic forces as well as the mass forces and their 
distribution from the calculation. But these forces are 
far from negligible as far as the magnitude of ph^se angle 
is concerned. To illustrate: It can be proved that with 
two degrees of freedom - wing bending and wing torsion - 
flutter is possible only when the product 

m s > (20) 
J. o 

whereby m is the wing mass per unit length of the span 
and s the distance of the e.g. of the wing element behind 
the quarter-chord point. Applied to a cantilever wing 
without aileron, this simply means shifting the e.g. of the 
individual wing sections near the quarter-chord points, in 
order to prevent flutter at any air density p # Such wings, 
although with ailerons, are found on the M 20, . 

The majority of flutter cases described In -sept ion III, 
probably could have been prevented by careful mass balanc- 
ing of the ailerons. This method, originally pointed out 
by von Baumhauer. and Soiling (reference 3) in 1923, has fre- 
quently been discussed since then in the literature. It is 
a fact, however, that practically all older models had un- 
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"balanced ailerons and that flutter was less frequent and le 
less dangerous at the then comparatively low flying speeds* 

On the other hand, mass-balancing the ailerons after 
an airplane is "built, means an expense of weight, requir- 
ing up to 0.5 percent of the airplane weight aside from 
increased drag. Quite obviously, subsequent modification 
of all existing types was therefore out of the question, 
particularly as the need for such measures did not seem 
very apparent as far as ,the older types were concerned. 

For the new types the DLA at first recommended mass 
balancing; later this. was incorporated into the airplane- 
design specifications. An aileron originally designed for 
mass balance is not much heavier than one designed without 
balance. The purpose of mass balancing is to reduce the 
aileron amplitude to harmless magnitude in all modes that 
may cause flutter. For modes such as shown in figures 22 
and 23, even a "partial mass balance" may be all that is 
required. 

In other cases, however, it is necessary to effect a 
complete and careful balance because it requires a ten- 
times-greater aileron amplitude in order to get out. of the 
range of the minimum of critical speed, according to the 
diagrams. This fact has not always been sufficiently rec- 
ognized. 

The success of these preventative measures should be 
checked on the oscillation bench, because even a complete- 
ly balanced aileron may oscillate. Possible causes are: 

1. Mass coupling due to the co-oscillating air mass, 

particularly when the gap between wing and aile- 
ron is small and the aileron is not aerodynamic- 
ally balanced, 

2. Kinematic coupling with complicated and indistinct 

static structure of the cellule. 

3. Lack of torsional stiffness of ailerons. 

4. Natural oscillations of the system: left aileron, 

controls, right aileron. 

The latter oscillation is particularly dangerous 
when coincident with a symmetrical natural mode of the 
wings. In the vicinity of the resonance point, the phase 
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angle of the aileron notion relative to the wing motion 
changes profoundly with the frequency, so that the phase 
angle most favorable for flutter can easily occur. Reme- 
dies are: changes in the natural frequency of the aileron 
control or artificial damping of the aileron deflection. 

Some of the above cases prove, at any rate, that even 
ailerons nearly or completely mass- balanced may occasion- 
ally develop flutter, in which case, however, the reduced 
frequency seems to he "below the average, i.e. 0) < 0.9. 

If the aileron were to "be considered as nonexisting 
so far as flutter is concerned, they would have to have 
not only zero amplitude on the oscillation "bench hut also 
complete aerodynamic "balance - at least, within the range 
of small angles of attack and aileron deflection. It can 
he proved theoretically that otherwise the circulation may 
cause a purely aerodynamic coupling, which lowers the crit- 
ical speed relative to that of the wing flutter in hending 
and torsion. 

If all these conditions for preventing aileron oscilla- 
tion were fulfilled, the displacement of the e.g. axis of 
the wing near to the quarter-chord point, would practically 
suffice for flutter prevention. Obviously, this is predicat- 
ed on the assumption that the two-dimensional theory of wing 
flutter is substantially correct, which cannot be summarily 
taken for granted with complicated wing shapes. The wing 
stiffness of such an airplane could be arbitrarily low, 
provided no other lower stiffness limits existed. 

Such flutter prevention, however, requires a large 
number of design changes of such a radical nature that in 
many cases it would be tantamount to a new departure in 
design methods. In view of this fact, it seemed more ex- 
pedient to increase the wing stiffness as long as consist- 
ent with minimum weight. This was the reason why this pre- 
ventative measure was resorted to at first. Greater stiff- 
ness leads to higher wing frequencies and consequently 
higher critical speed. As a result, the flutter, while not 
altogether prevented is, however, moved up into a speed 
range above the highest speed which can be reached, A con- 
tributing factor was the consideration that the wing it- 
self must have a certain minimum stiffness in order to pre- 
vent static torsional instability and reversal of aileron 
effect at high flying speed. 
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With further increase of speed, however, a point is 
reached where the simple expedient of increased stiffness 
is no longer compatible with the weight and where it will 
"be necessary to combine all known measures for the pre- 
vention of flutter. 



Translation by the national Advisory 
Committee for Aeronautics. 
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Figs. 1,2,3,4,5,6 




with servo-rudder Figure 2.- Wing section 



.Strut 




Figure 3.- Critical fonn of 

oscillation of the 
He 8 at 550/m in. scale of 
wing contour 1:60 




Figure 4.- Dangerous oscillation 
of L 78 at 860/min, 
scale 1:80 




Lower wing 
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Figure 5.- Dangerous oscillation 
of He 60 780/min, 
scale 1:80 



Figure 6 # - Dangerous oscillation 
of He 46c at 815 /min, 
scale 1:60 
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Tigs. 7,8,9,10,11,12 



N Strut 




Figure 7.- Dangerous oscillation 
of KL 1 A at 675/min. 
scale 3:200 



Figure 8.- Dangerous oscillation 
of Do lOj scale 1:100 

Lines of break on left and 

right wing 
— j — Break or stretch of cahles 




Figure 10.- Dangerous oscillation 

of AC 12 E; scale 3:200 



Figure 9.- Dangerous oscillation 
of L 102 at 835/min* 
scale 1:80 
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Figure 11.- Dangerous oscillation 
of Do 12 at 580/min 
asymmetrical excitation 
scale 3:100 




Figure 12.- Dangerous oscillation 
of U 28 at 770/min. 
asymmetrical excitation 
scale 1:100 
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*igs. 13,14,15,16,17,18 




Lower wing 

Figure 13*- Dangerous oscillation 

of S 24 at 1215/min. 
scale 1:150 




Figure 14,- Width 
of res- 
onance 
curve 
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Figure 16.- Aileron 
vibration 




Figure 15*- Plate 
vibration 
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Figs. 19.20,21, 22,23 





Figure 19.- Wing torsion 
and aileron 
motion, 
r * 0.25 



Figure 21.- Wing torsion 

about l/4 
chord point and aileron 
motion, t * 0.25 



Figure 22.- Wing torsion 

about 3/4 
chord point and aileron 0 
motion, t * 0.15 
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Figure 20.- Wing torsion 

about l/4 
chord point and aileron 
motion, t * 0.15 



Figure 23,- Wing torsion 

about 3/4 
chord point and aileron 
motion, t • 0.25 
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